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Reconstruction of density matrices is important in NMR quantum computing. An analysis is 
made for a 2-qubit system by using the error matrix method. It is found that the state tomography 
method determines well the parameters that are necessary for reconstructing the density matrix in 
NMR quantum computations. Analysis is also made for a simplified state tomography procedure 
that uses fewer read-outs. The result of this analysis with the error matrix method demonstrates 
that a satisfactory accuracy in density matrix reconstruction can be achieved even in a measurement 
C^l ' with the number of read-outs being largely reduced. 

o 

Q 

The study of quantum computers has attracted considerable attention since Shor in 1994 introduced a quantum 
mechanical algorithm for efficient factoring of large numbers jjj . In another remarkable work, Grover in 1996 discovered 
s ! \ that quantum mechanics can help to speed up data search in an unsorted database |^J^]- Quantum mechanics that, 
for nearly a centenary, has been a basic tool for understanding the microscopic world is becoming a powerful new 
^*T* ' weapon for computation, communication and information-processing. 

How to realize quantum computing experimentally has sparked an explosion of interest. Among many proposed 
' physical systems to implement quantum computation, such as trapped ions, optical photons, quantum dots, and so 
y—i , on, the NMR quantum computation is particularly attractive because nuclear spins are extremely well isolated from 
■ their environment and readily manipulated with modern NMR techniques. Comprehensive algorithm realizations 
have been accomplished for the Deutsch-Jozsa algorithm in 2-qubit and 3-qubit systems j§-[3| and a 5-qubit system 
0, and for Grover's algorithm with 2-qubit ||[| and 3-qubit systems [|l0). The order-finding problem in a 5-qubit 
O |' system has recently been demonstrated in NMR |ll|] , and a Cat state is prepared in a 7-qubit system jl2] . 

The fundamental elements for information-processing in NMR are two-level nuclear spins that are bound together 
in a single molecule. Put forward in 1997 by Cory et al. |l3| and Gershenfeld et al. JO) , the NMR scheme uses 
bulk numbers of molecules. The NMR techniques cannot control the quantum states of individual molecules; instead, 
all the molecules in the sample are manipulated in parallel. In fact, a liquid NMR sample is initially in a thermal 
equilibrium at room temperature. The directions of the nuclear spins have a Boltzman distribution and are not 
polarized along the strong magnetic field. 

Reconstruction of quantum state, or quantum state tomography, is not only important in quantum computation, 
but also important for general purposes. It has attracted many investigations. For instance, discrete Wigner function 
has been used to infer the quantum states of finite-dimensional systems from measurement by Leonhardt ]l5| , ^6[ . 
Walser et al have proposed a scheme to reconstruct the state of a single quantized cavity mode ||l7|] . Amiet and 
Weigert have studied various aspects of reconstruction of density matrix for a spin state in pure or mixed states 
JTs] . Reconstruction of state of light has been extensively studied by Leonhardt |21[ . 

In NMR quantum computation, the liquid ensemble is described by the density matrix. The state of the NMR 
computer can be obtained by state tomography technique |22]| . In order to extract the density matrix, for example, 
for a 2-qubit system, 18 read-outs have to be performed. In general, for an n-qubit system, construction of the density 
matrix requires 3™xn read-outs. After signal read-out, the area of the spectrum is integrated and the density matrix 
is reconstructed through numerical methods. Obviously, the amount of work in experiment and post-processing is 
huge when n becomes moderately large. 

In this paper, we study the reconstruction of density matrix in NMR. The mathematical structures of the density 
matrix reconstruction is analyzed. We found that the full state tomography read-outs, which consists of 18 read-outs 
for a 2-qubit system, is well conditioned. However, a smaller number of read-outs can also determine the density 
matrix well. We compared the results of 2-qubit density matrix in NMR experiment with the full read-outs and a 
much reduced read-outs respectively and found that they agree with the theoretical analysis. 

The paper is arranged as follows. In Section II, we describe the procedures to construct the density matrix. We 
analyze the density matrix reconstruction procedure in section III. In Section IV, we discuss the reduction of number 
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of read-outs in constructing the density matrix. We apply the analysis to a 2-qubit NMR experiment and the results 
indicate that it is possible to construct a density matrix with much reduced number of read-outs and at the same the 
loss of accuracy is small. Finally in Section V, we give a summary. 



II. RECONSTRUCTION OF DENSITY MATRIX 



In an NMR measurement, each read-out pulse can only give some off-diagonal matrix elements of the density matrix. 
To obtain the rest matrix elements, one has to rotate the original density matrix through rotational operations. In 
a 2-qubit system, in order to construct the density matrix, one needs to perform the following operations P2f : II, 
IX, IY, XI, XX, XY, YI, YX, and YY. Here, I, X and Y stand for, respectively, the identity operation, a 90 degree 
rotation about the x-axis, and a 90 degree rotation about the y-axis. Thus, in a state tomography, these operations 
are performed before NMR read-out measurements. We explain this through an example below. 

Suppose that we use the nuclear spins of H and P in a phosphorous acid as our qubits. Since for a usual NMR 
system, only one nuclear spin can be measured at a time, we have to perform the measurement separately for the 
two nuclear spins. We first start a computation and do a measurement on H at the required stage. We then restart 
the computation from the beginning and measure the signal corresponding to P. Next, we restart the computation, 
but this time perform the operation IX at the required stage before measuring the signal. This process is separately 
carried out for H and P nuclear spins, and the nine operations are successively performed for each of them. This 
means that we have to perform totally 9 x 2 = 18 read-outs. 



Taking |) = |1) and | J,) = |0), the matrices of the nine operations are: 
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Let us assume that the density matrix of the system takes the form 
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(i) 



Pll Pl2 Pl3 Pli 

Pl2 P22 P23 P2i 

Pl3 P23 P33 P3A 

Pli P2A Pli P44, 



Xi 2 2 + IXll %3 + iXl2 %i + «£l3 

X2 — ixn 25 x§ + ixi4 x-j + ixis 

X3 - 1x12 x 6 - ixu x 8 x 9 + ixie 

24 — 1x13 X-; — ixis x 9 — xio 



(2) 



The NMR read-out signal can only give X2 + ix\\ and xg + ixie for the nuclear spin of P, and 23 + ixu and 27 + 1215 
for the nuclear spin of H. For example, the element P13 in Eq. (0) corresponds to the left peak of the spectrum of H, 
while the element P24 to the right peak of it. Similarly, the element pi2 corresponds to the left peak of the spectrum 
of P, and the element P34 to the right peak of the same spectrum. To obtain other elements in the density matrix, we 
have to perform one of the nine operations for the system so that the desired elements are transformed to the positions 
labeled as 12, 13, 24 and 34 in the density matrix, and thus can be measured. The read-out gives the elements p' 12 , 
p' 13 and p'24 , P34 in the rotated density matrix. These rotated matrix elements are linear combinations of the original 
matrix elements. For each of the nine operations, one makes two measurements for the nuclear spins of H and P, and 
each measurement provides two matrix elements which contains a real and an imaginary part. Altogether, we finally 
obtain 4 x 9 x 2 = 72 equations with 16 unknowns. 

The coefficients in these equations form a matrix with the size of 72 x 16. The rank of this coefficient matrix is 
only 15 since the matrix can be added with a matrix of a constant times the unit matrix without changing the results. 
In practice, we can add one more equation by letting the trace of the density matrix be 1 after normalizing the 
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integrations of the spectrum. Thus, we have totally 73 equations. To determine the elements in the density matrix, 
one needs to solve the following set of linear equations for xc 

16 

^A ai x i = B a , {a =1,2,... 73}, (3) 

i=l 

where A a i is the coefficient of Xi in the a-th equation and it varies with read-out for different rotations, and B a is the 
integrated area of the spectrum. There are certainly redundant expressions in (^|) since the number of equations is 
more than the number of unknowns. The standard way of dealing with this problem is to use the least square fitting 
procedure that is widely used in various problems in science and engineering. We minimize the quantity x 2 defined 
as 

16 

x 2 = E(E^^-^) 2 - ( 4 ) 

a i=l 

To find the minimum, we carry out a variation procedure on \ 2 with respect to all parameters, which gives 

16 

yijXj — bi, (5) 



E^ 



3=1 

where 

Cjj = A ai A aj , and bi = ^B a A al . (6) 

a a 

The number of equations in (Q) is now equal to the number of unknowns. These linear equations can be solved 
by standard numerical method such as the Gaussian elimination. In principle, such problems can be solved in this 
way. However, in most of our cases, not every parameter can be well determined in experiment. In the least square 
fitting procedure, these less well-determined matrix elements might possibly bring spurious values into the numerical 
calculation while leaving the x 2 value small. In other words, the x 2 might not be very sensitive in response to a 
big variation in experiment. For instance, if we have a set of equations, x+y=cl, x+1.001y=e2, this set of equations 
determines the sum of x and y very well, and the difference between x and y is poorly determined. The x 2 is insensitive 
to a change in x-y while leaving x+y unchanged. 



III. AN ANALYSIS OF THE DENSITY MATRIX RECONSTRUCTION 

In order to obtain the density matrix, one needs to use the least-square fitting method [p3| to derive them from 
the experimental data. Similar problem occurs actually in other fields of science as well. Here, we adopt a variant 
of the least square-fitting method from nuclear physics to solve the problem. In the large-scale nuclear shell model 
calculations, Wildenthal et al. employed the error matrix method pl| ] to analyze sensitivity of the nuclear structure 
to experimental input. In their method, instead of solving Eq. (|^) directly, one first solves the eigenvalue problem for 
C defined in (|), 

UCU+ = C d . (7) 
Here, U is the unitary matrix that diagonahzes C. Cd is called the error matrix (24j. Then Eq. (||) becomes 

C d y = b', (8) 
where yi — UijXj and b i = Uijbj. b { contains experimental information. Since Cd is diagonal, we can determine j/j by 

Vi = Tf^T- ( 9 ) 

[ydjii 

b i contains experimental uncertainties, and any change in b i will cause yi to change. However, if the diagonal matrix 
element {C d )u is large (say, of the order of 1), yi will be insensitive to changes in b it and thus it can be well determined. 
Conversely, a small {Cd)u (say, 0.001) means that the corresponding yi is very sensitive to experiment, and any small 
variation in b' will cause a big change in y. In this case, y is not well determined, and special effort is necessary to 
ensure that b' is sufficiently accurate. If this is not possible, then during the fitting process, the y's that corresponds 
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to small eigenvalues of C will be kept constant by physical considerations. In nuclear structure studies, the critical 
constant value was chosen as 0.001. 

The C matrix for a 2-qubit system with all the nine transformations performed (with 18 read-outs, and 73 linear 
equations) is 
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Solving the eigenvalue problem, and we obtain the 16 eigenvalues of C. They are: 4, 4, 4, 4, 4, 6, 6, 3, 3, 2, 6, 4, 4, 
4, 4, 6. The j/'s are the combinations of the x coefficients 

Vl = -0.37a; 3 + 0.82a;4-0.24a; 6 ; 

?/2 = -0.19xi + 0.41x2 - 0.17x 3 - 0.32x 4 - 0.19x 5 - 0.59x 6 + 0.17x 7 - 0.19x 8 - 0.41x 9 - 0.19xi ; 

y 3 = -0.0087xi - 0.093x2 - 0.57x 3 - 0.39x 4 - 0.0087x 5 + 0.42x 6 + 0.57x 7 - 0.0087x 8 + 0.093x 9 - 0.0087xi ; 

?/4 = 0.30xi - 0.30x2 - 0.093x3 - 0.26x 4 + 0.30x 5 - 0.61x 6 + 0.093x 7 + 0.30x 8 + 0.30x 9 + 0.30xi ; 

2/5 = -0.35xi - 0.48x2 + 0.027x 3 - 0.035x 4 - 0.35x 5 - 0.20x 6 - 0.027x 7 - 0.35x 8 + 0.48x 9 - 0.35xi ; 

2/6 = -0.55x2 + 0.44x 3 + 0.44x 7 - 0.55x 9 ; 

2/ 7 = -0.44x 2 - 0.55x 3 - 0.55x 7 - 0.44x 9 ; 

2/8 = 0.71x5 - 0.71x 8 ; 

2/ 9 = -0.71.xi + 0.71xi ; 

2/io = 0.50xi - 0.50x 5 - 0.50x 8 + 0.50xi ; 

2/11 = 0.71xi2 + 0.71xi 5 ; 

2/i2 = -0.71xl2 + 0.71xi5; 

Vi3 = x 13 ; 

2/14 = X14,; 

2/i5 = OJlxii - 0.71xi 6 ; 

2/i6 = 0.71xii + 0.71xi6. 

Some of the y's are directly the x parameters in the density matrix such as 2/13 and 2/14, and others are the combinations 
of the x parameters. They are all well determined by the experimental data. 



IV. REDUCTION OF EXPERIMENTAL READ-OUTS 



In the process of the density matrix reconstruction discussed above, each signal read-out provides us two equations. 
By performing all the nine operations for P and H, and plus the normalization condition, we have 73 equations. 
However, these equations are over-determined. The question is if it is possible to determine the density matrix with 
less read-outs. It will be very interesting to see what is the minimum number of operations to determine the density 
matrix without loss of much accuracy. Intuitively, we may think about 4 read-outs because there are altogether 16 
unknowns. However, a detailed analysis of the rank of equations indicated that any 4 read-outs combined together 
can not provide sufficient numbers of independent equations. The minimum number of read-outs to be performed is 
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5. In table (|) we give the various set of 5 read-outs. In the table, operation 1 is the II with H signal acquisition(see 
table caption for details). From the table we see that not arbitrary combination of 5 operations are possible. 

To see how well the density matrix is determined by a smaller number of read-outs, we choose a 6 read-outs case. 
The six read-outs are: II, IX, IY, XX for one of the nuclear spins (H), and II, IX for the other (P). In this situation, 
the corresponding C matrix becomes: 
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The eigenvalues of C matrix are: 1, 2, 1, i, i, 1, 4, 2, 4, i, 1, 3, 1, 2, 1, 2. It can be seen that the eigenvalues are all 

quite big. The y expressions corresponding to the eigenvalues are: 

2/i = -X6l 

2/2 = x 2 ; 

V3 = x 4 ; 

Vi = -0.79x1 + 0.21x5 + 0.58x8; 

2/5 = -0.21xi + 0.79x 5 - 0.58x 8 ; 

y<3 = 0.71x 3 - 0.71x 7 ; 

2/7 = 0.71x 3 + 0.71x 7 ; 

Us = x 9 ; 

2/9 = -0.50xi - 0.50x5 - 0.50x8 - 0.50xi ; 

2/io = -0.29xi - 0.29x5 - 0.29x8 + 0.87xi ; 

2/11 = -0.71xi2 + 0.71xi 5 ; 

2/12 = 0.71xi2 + 0.71xi 5 ; 

2/i3 = -0.71xi3 - 0.71xi 4 ; 

2/14 = -0.71xi3 + 0.71xi 4 ; 

2/i5 = -0.71xii - 0.71xi 6 ; 

2/16 = -0.71xn + 0.71xi 6 . 



(12) 



The density matrix elements are as well determined as the one with 18 read-outs. However, the saving in the numbers 
of read-outs is great (18 — 6 = 12). Because we have less read-outs now, statistics in this case is surely poorer, and this 
will reduce accuracy of the parameters determined. Nevertheless, because of the intrinsic mathematical structure, the 
uncertainties in the experimental data affect the parameters insensitively, and the loss in accuracy is not very big. 

As an example, we have analyzed the density matrix for a 2-qubit system in an NMR experiment P6[ . The 
theoretical prediction of the density matrix is 



Pth = 



-0.063 



0.31 
0.31 
0.31 
- 0.13i 



0.31 
0.31 
0.31 

-0.063 + 0.13i 



0.31 
0.31 
0.31 

-0.063 + 0.13i 



-0.063-0.13* 
-0.063 -o.m 
-0.063 -0.13i 
0.063 



(13) 
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When constructing the density matrix using all 73 equations, we obtain 



Pall 



0.36 0.33 - 
0.33 + 0.087* 

0.31 - 0.037* 0.28 + 

-0.034 + 0.22i -0.079 



0.087* 0.31 + 0.037* -0.034 - 0.22* 

0.30 0.28 - 0.022* -0.079 -O.Ui 

0.022* 0.23 -0.044 -0.13* 

h0.14* -0.044 + 0.13* 0.12 



(14) 



If we take only 49 equations out of the 73 equations, which means that we take only 12 read-outs (which contains all 
9 read-outs for H and the II, IX, IY read-outs for P) instead of the complete 18 read-outs, the density matrix is 



0.31 
0.29- 
-0.074 



0.37 
f 0.13* 
0.034* 
+ 0.17* 



0.31 - 0.13* 
0.28 

0.37-0.0091* 
-0.059 + 0.Ui 
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-0.074 
-0.059 
-0.058 



-0.17* 
- 0.14* 

-o.m 

0.094 



(15) 



We see that the experimental density matrices arc nearly identical for the two cases. Defining the error measurement 



5 = 



iPexp — Pth||2 



(16) 



where |||| is the norm of a matrix, we get the error of p a n relative to pth 17%, and the same error of p±2 to pth, which 
roughly says that p a ii or p\2 is about 83% of pth- Thus we have seen that the number of read-outs is reduced quite a 
lot, but the accuracy in the density matrix is not much influenced. For the extreme case with 6 read-outs only, the 
density matrix is: 



Pa 



0.52 0.29-0.21* 0.29 + 0.068* 

0.29 + 0.21* 0.16 0.37 + 0.085? 

0.29- 0.068* 0.37- 0.085* 0.32 

0.025 + O.Ui -0.045 + 0.13* 0.0039 + 0.15* 



-0.025 -O.Ui 
-0.045 -0.13* 
0.0039 - 0.15* 
0.18 



(17) 



Still, it is quite close to pth- The relative error of pg to pth is 32%. We see that, although the number of equations 
is reduced more than a half, the errors do not increase as much as one would think. Since the reconstruction of the 
density matrix in NMR quantum computing is so tedious, the present work suggests a way to reduce the number of 
read-outs if the accuracy is not so highly required. 



V. SUMMARY 



In conclusion, we have analyzed the error matrix for the density matrix reconstruction in the 2-qubit NMR quantum 
computing. We have found that the number of read-outs can be reduced greatly without significant loss in the accuracy. 
Our analysis can be easily extended to NMR systems with a larger number of qubits. 
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the Fok Ying Tung Science Foundation. 
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